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Phase separation in a two species Bose mixture 
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We obtain the ground state quantum phase diagram for a two species Bose mixture in a one- 
dimensional optical lattice using the finite size density matrix renormalization group(FSDMRG) 
method. We discuss our results for different combinations of inter and intra species interaction 
strengths with commensurate and incommensurate fillings of the bosons. The phases we have 
obtained are superfluid, Mott insulator and a novel phase separation, where the two different species 
reside in spatially separate regions. The spatially separated phase is further classified into phase 
separated superfluid(PS-SF) and Mott insulator(PS-MI). The phase separation appears for all the 
fillings we have considered; whenever the inter-species interaction is slightly larger than the intra- 
species interactions. 
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I. INTRODUCTION 

Studies of quantum phase transitions are currently of 
great interest as they provide important insights into a 
wide variety of many-body systems [l], Q ■ The pioneer- 
ing observation of the superfluid (SF) to Mott insulator 
(MI) transition in an optical lattice using cold bosonic 
atoms [U, which had been predicted by Jaksch et. al. 
[H, highlights the exquisite control of the inter atomic 
interactions that is possible in such systems. In that ex- 
periment, performed using 87 i?6 atoms, the tunneling of 
the atoms to neighboring sites and also the strength of 
the on-site interactions were controlled by tuning and/or 
detuning the laser intensity in order to achieve the transi- 
tion from the SF phase (random distribution of atoms) to 
the MI phase where there are a fixed number of atoms per 
site Q . Recent developments involving the manipulation 
of ultracold atoms have led to the realization of genuine 
one dimensional systems such as the Tonks- Girardeau 
gas Q ■ Several interesting phenomena including the SF- 
MI transition have been observed in one-dimensional op- 
tical lattices [f|. 

In the past few years, on the theoretical side, many in- 
vestigations have been carried out using a single species 
of bosonic atoms in optical lattices @, [1[ . Recently cold 
bosonic mixtures [9], fermions 

.10] 

and Bose-Fermi mix- 
tures [Til [l2j in optical lattices have attracted much at- 
tention. Mixtures of different species are very interesting 
since additional phases could appear due to the inter- 
species interactions [13, LLJ] • 

In the present work, we consider a system with two 
species of bosonic atoms or equivalently, bosonic atoms 
with two relevant internal states. The two species shall be 
called a and b type respectively. The low-energy Hamil- 



tonian is then given by the Bose-Hubbard model for the 
two boson species: 



H = - ta E ^ + h -°) - tb E + h -°) 

<ij> <i,j> 

+y a E<«- 1 ) + ?E**- 1 ) 



+u ab j2 n ?? 



(i) 



*E-mail: tapan@iiap.res.in 
t E-mail: rvpai@unigoa.ac. in 
t E-mail: das@iiap.res.in 



Here cii (6j) is bosonic annihilation operator for bosonic 
atoms of a (b) type localized on site i. nf — ajflj and 
n\ = b\bi are the number operators. t a (t b ) and U a (U b ) 
are the hopping amplitudes between adjacent sites (ij) 
and the on-site intra-species repulsive energies, respec- 
tively for a (b) type of atom. The inter-species interac- 
tion is given by U ab . The hopping amplitudes (t a , t b ) 
and interaction parameters (U a , U b , U ab ) are related to 
depth of optical potential, recoil energy and the scatter- 
ing lengths. The ratio JJ ab /U a ' b can be controlled to a 
wide range of values (l5j experimentally. In this work 
we consider inter-exchange symmetry a <-> 6, implying 
t a = t b = t and U a = U b = U and study the effect of 
inter-species interaction on the ground state of model([T]) 
in one-dimension. We set our energy scale by t = 1. 

The model ([1]) has been studied earlier using the mean 
field [Tp , Monte-Carlo [l3[ and the Bosonization meth- 
ods [14| and this has resulted in the prediction of the 
basic structure of its ground state phase diagram. The 
Bosonization study predicts phase separation (PS) for 
large values of the inter-species interaction U ab by con- 
sidering one species of bosons to be hard core and the 
other to be in the intermediate to hard core regime [l4l |. 
Phase separation has also been found using a variational 
method based on the multi orbital best mean field ansatz 
[l6| . However, a clear picture of the transitions pertain- 
ing to the SF, MI and PS phases has not emerged so far. 
In order to achieve this, we consider the influence of the 
inter-species interaction U ab on these phases, by carry- 
ing out a systematic study of its effect on the ground 
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state of model (QJ in one-dimension using the finite size 
density matrix renormalization method [H, [13 • The fol- 
lowing part of the paper is organized as follows. Section 
II contains the details of our finite size density matrix 
renormalization method. Section III contains our results. 
We end with concluding remarks in Section IV. 



II. FSDMRG CALCULATION 

The Finite-Size Density-Matrix Renormalization 
Group (FSDMRG) method has proven to be very useful 
in studies of one-dimensional quantum systems 0, [l7| . 
The details of this method are given in a recent review 
by Schollwock [18[ . The open boundary condition is pre- 
ferred over periodic boundary condition for this method 
because the loss of accuracy which increases with the size 
of the system is much less in the former than the latter. 
In the conventional FSDMRG method, the lattice is first 
built to the desired length (L) using the infinite system 
density- matrix renormalization method (DMRG). The 
finite size sweeping is done only for this desired lattice 
size L. We use a slightly modified form of the FSDMRG, 
where we sweep at every step of the procedure and not 
just for the case which corresponds to the largest value 
of L. This enables one to obtain accurate correlation 
functions. Furthermore, since the superfluid phase in 
models such as Eq. {TJ, in d = 1 and at T = 0, is 
critical and has a correlation length that diverges with 
the system size L, finite-size effects must be eliminated 
by using finite-size scaling as we show later. For this 
purpose, the energies and the correlation functions, 
obtained from a DMRG calculation, should converge 
satisfactorily for each system size L. It is important, 
therefore, that we use the FSDMRG method as opposed 
to the infinite-system DMRG method especially in the 
vicinities of continuous phase transitions. 

In the FSDMRG method the bases used for left- and 
right-block Hamiltonians are truncated by neglecting the 
eigenstates of the density matrix corresponding to small 
eigenvalues which leads to truncation errors. If we re- 
tain M states, the density-matrix weight of the discarded 
states is Pm = Sa=i(l — w a), where u> a are the eigenval- 
ues of density matrix. Pm provides a convenient measure 
of the truncation errors. We find that these errors depend 
on the order-parameter and correlation length for a given 
phase. For a fixed M, we find very small truncation er- 
rors in the gapped phase and the truncation errors are 
largest for the SF phase. In our calculations we choose M 
such that the truncation error is always less than 5 x 10 -5 
and we find that M = 128 suffices. 

The number of possible states per site in the model([T]) 
is infinite since there can be any number of a and b species 
bosons on a site. In a practical FSDMRG calculation 
we must truncate the number of states n max allowed per 
site. The value of n max , of course, will depend on the on- 
site interaction U. The smaller the value of U the larger 
must be n max . From our earlier calculation [f| on related 



models, we find that n max = 4 is sufficient for the value 
of U considered here. This implies, for modelfl]), 4 states 
each per site for a and b species bosons and a total of 16 
states per site. This corresponds to a truncation of bases 
of left (right) block from 16M to M in each FSDMRG 
iteration. 

Before proceeding further we give a brief summary of 
our results. The various parameters that we calculate to 
study the ground state properties of model ([1]) are the 
energy gap Gl , which is the difference between the ener- 
gies needed to add and remove one atom from a system 
of atoms, i.e., 

G L = E L (N a + 1, N b ) + E L (N a - 1, N b ) - 2E L (N a , N b ) 
and the on-site density correlation function 

( n f) = (lpOLN a N b \nf\lp LN a N b )- (3) 

Here a, is an index representing type a or b bosons, 
EL{N a , Nb) is the ground-state energy for a system of 
size L with N a (N b ) number of a (b) type bosons and 
\i>0LNa,N h ) is the corresponding ground-state wavefunc- 
tion, which arc obtained by the FSDMRG method. In 
d = 1 , the appearance of the MI phase is signaled by the 
opening up of the gap Gl^oc- However, Gl is finite for 
finite systems and we must extrapolate to the L — > oo 
limit, which is best done by using finite-size scaling [f|. 
In the critical region, i.e., SF region, the gap 

GL^L-ifiL/S), (4) 

where the scaling function f(x) ~ x, x — > and £ is 
the correlation length. £ — > oo in the SF region. Thus 
plots of LGl versus U, for different system sizes L, con- 
sist of curves that intersect at the critical point at which 
the correlation length for L = oo diverges and gap Goo 
vanishes. 

Defining the ratio of the inter and intra species in- 
teractions A = U ab /U, we study the ground state of 
model flU for A < 1, A = 1 and A > 1. The ground 
state exhibits some similarities as well as differences in 
each of the cases. When the kinetic energy is the domi- 
nant term in the model, the ground state is in 2SF (both 
a and b species are in the SF phase) state for all A. This 
similarity is, however, lost when the interactions domi- 
nate. For A < 1, i.e., U ab < U, the large U phase is 
Mott insulator with non-zero energy gap in the ground 
state. This state has an uniform local density of bosons 
for each species, i.e., = (n b ) for all i. The 2SF 

to MI transition is possible only when the total density 
p = p a + Pb is an integer. For U ab - U, the 2SF-MI tran- 
sition for model ([T|) is similar to the SF-MI transition for 
single species bosons with the same density of bosons. 
For A > 1 and for small values of U, the ground state 
is a 2SF state. However, when U increases, the ground 
state first goes into superfluid phase with a and b bosons 
spatially separated into different regions of the lattice. 
This is the case when p a — p b — 1/2. This phase may 
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be called the phase separated superfluid(PS-SF). There 
is no gap in the ground state energy spectrum and the 
phase separation order parameter defined as 



O 



PS 



(^0LN a N b \(K ~ r$\)\ll>0LN a N h )- (5) 



When the densities 
p b = 1/2, no PS-MI 



is non-zero. A further increase in U results in opening up 
of the gap in the energy spectrum. This Mott insulator 
has a non-zero phase separation order parameter and it 
may be called the phase separated Mott-Insulator(PS- 
MI). The total local density (rii)(= ((nf + n b ))) = p 
remain uniform across the lattice 
are different, for example p a = 1. 
is found and the ground state has only 2SF and PS-SF 
phases. When p a = 1, pt, = 1 we find, for A = 1.05, no 
PS-SF phase and the transition is directly from 2SF to 
PS-MI. We now present the details of our results. 



III. RESULTS AND DISCUSSIONS 

In the absence of the inter-species interaction U ab , 
model ([I]) is an independent mixture of the individual 
species of bosons. So the nature of the ground state of 
model JT]) depends only on the density of the individual 
species of bosons: p a , pb and U a = U b — U , the on-site 
interactions. For example, if p a =/= n, pt =/= n, where n is 
an integer, the ground state is always in the superfluid 
phase irrespective of the strength of the on-site interac- 
tion U. The Mott insulator is possible only when either 
p a = n or pb = n. Based on the values of p a , pb and U, 
the ground state is categorized as 2SF(both a and b type 
bosons in SF phase), SF+MI(a boson in SF and b in MI 
phase or vice-versa) and MI+MI(both a and b bosons in 
MI phase). For p a ^ n, pb + 1 n or for any values of p a , 
Pb, but U < U Cl where U c is the critical on-site interac- 
tion for SF to MI transition, the ground state is always 
in the 2SF phase. SF+MI phase is possible for p a ^ n, 
Pb = n ( or vice- versa) and U > Uc ■ If both p a = Pt = n 
and U > Uc, we have MI+MI phase. In order to in- 
vestigate the influence of U ab on these ground states, we 
consider three cases: A < 1, A = 1 and A > 1, where 
A = U ab /U. In each of these three cases, we consider 
three different ranges of densities:(i) p a = Pb = 1/2, (ii) 
p a = l.pb = 1/2, (iii) p a = pb = 1. The choice of these 
three cases are made to understand the effect of the inter- 
species interaction on 2SF, SF+MI and MI+MI phases. 
We now discuss each case below. 

(i)p a = Pb = 1/2: 
As discussed in the previous paragraph, for this case, 
there is no MI phase if U ab = and the model (JI) will 
have only the 2SF phase. However, with the introduction 
of inter-species interaction, the 2SF phase is destroyed. 
For example, Figure |T]) shows a plot of scaling of gap 
LGl versus U for A = 1. Curves for different values of L 
coalesce for U < U c — 3.4 indicating a gapped MI phase 
for U > U c . The emergence of this phase is due to the 



intra-species as well as interspecies interaction strengths. 
The fact that U c — 3.4, indicates that the model (JTJ) 
when A — 1 behaves like a single species of bosons at 
unit density @. These results are along the expected 
lines because, when U ab — U, every boson in the system 
interacts with rest of bosons, irrespective of whether they 
are of type a or b, with the same strength and therefore 
the species index become irrelevant. However, the situa- 
tion changes when the inter-species interaction U ab ^ U. 
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FIG. 1: Scaling of gap LGl is plotted as a function of U 
for different system sizes for p a = pt, = 1/2, A = 1. The 
coalescence of different curve for U < 3.4 shows a Kosterlitz- 
Thouless-type 2SF-MI transition. This transition is similar to 
the SF-MI transition for single species Bose-Hubbard model 
for p= 1 [|. 

For A < 1, i.e, U ab < U, the system still undergoes 
2SF-MI transition when the on-site repulsion increases, 
but with a higher U c . For example figure ([2]) shows a 
plot of scaling of gap LGl versus U for A = 0.5. The 
critical U C (A = 0.5) ~ 5.4 is substantially greater than 
U C (A = 1) ~ 3.4. The ground state of the model © 
for p a = Pb = 1/2, A < 1 consists only of 2SF and MI 
phases. The transition from 2SF to MI is of Kosterlitz- 
Thouless-type. 

When A > 1, the scenario is drastically different from 
the one seen above. The on-site densities (nf) and (n b ) 
are plotted in Fig. for A = 1.05. It is clear from 
this figure that there is a spatial separation between the 
two different species of bosons for U = 4 and no spatial 
separation for U = 1. This highlights a phase separa- 
tion (PS) transition as a function of U. The question 
then arises whether this spatially separated phase is a 
superfluid or a Mott Insulator. In order to sort this out, 
we plot both the scaling of the gap LGl and the order 
parameter Ops for phase separation in Fig. (U|). It is 
evident from these figures that the transition to the MI 
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FIG. 2: Scaling of gap LGl is plotted as a function of U 
for different system sizes for p a = pt = 1/2, A = 0.5. The 
coalescence of different curve for U ~ 5.4 shows a Kosterlitz- 
Thouless-type 2SF-MI transition. 



FIG. 3: Plots of (n°) and (n-) versus i for U = 1 and U = 
4. These plots are for p a = Pb = 1/2, A = 1.05 and for 
system size L = 50. The deviation in (nf ) and (nj) near the 
boundaries for (7 = 1 is due to the open boundary condition 
used in our FSDMRG. 



phase happens at around U c ~ 3.4 and to the spatially 
separated phase around U c ~ 1.3. The gap remains zero 
for 1.3 < U < 3.4. Thus for the case p a — pb — 1/2 and 
A = 1.05, there are three phases: the superfluid phase 
(2SF) for U < 1.3, superfluid, but phase separated (PS- 
SF) for 1.3 < U < 3.4 and finally Mott Insulator, but 
again phase separated (PS-MI) for U > 3.4. It should 
be noted that the total density of bosons p = p a + Pb re- 
main constant through out the lattice, though bosons are 
space separated. The critical values of 2SF to PS-SF and 
PS-SF to PS-MI transition depends of the on the value 
of A. The detailed phase diagram in the A — U plane 
and the nature of the different phase transitions will be 
reported elsewhere. 

(n)p a = 1,^ = 1/2: 
In this case, when U ab = species a bosons undergo 
a superfluid to Mott insulator transition at U c ~ 3.4 
by virtue of having density p a = 1, while the b bosons, 
which has density p\, = 1/2, remains in the superfluid 
phase. However, when U ab < U, no transition from SF 
to MI was found for either of the two species of bosons. 
The Mott insulator phase of a bosons is completely lost. 
In the Fig. ((Ha)), we plot the length dependence of gap 
Gl for different U, which clearly indicates that the gap 
vanishes at L — > oo for all values of U considered. This 
emphasizes the fact that as far as the transition to the 
Mott insulator is concerned, when U ab < U, the total 
density must be an integer irrespective of the densities 
of the individual species of bosons; and it is this condi- 
tion that really matters. This condition remains same for 
jjab ^ jj Thus when the inter-species interaction is non 
zero as in the present case and the total density p ^ n, 
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FIG. 4: Plots of scaling of gap LGl (a) and order parameter 
for phase separation Ops (b) versus U demonstrate various 
phases for the case p a — pb = 1/2 and A = 1.05. 
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no Mott insulator phase is observed. 

The phase separation, however, happens when U ab > 
U. The local density distribution of different species of 
bosons are given in Fig. ©b)) for U = 1, 4 and A = 1.05. 
For {7 = 1, we find no phase separation, however, for 
C/ = 4, the a and b species bosons are phase separated. 
They rearrange in such a manner that the total density 
p = p a + Pb remains a constant. For example, when 
p a = 1 and p — 1/2, one-third of the region is occupied 
by the species b and two-third by the species a. The total 
density p being 3/2, the distributions of the a and b types 
of bosons follow the ratio of their densities. 
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This is consistent with the similar observations made for 
the case of p a = pb = 1/2. 
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FIG. 6: Scaling of gap LGl is plotted as a function of U for 
different system sizes for p a — pb = 1 and A = 1.0. The 
coalescence of different curve for U ~ 5.7 shows a Kosterlitz- 
Thouless-type 2SF-MI transition. 

The phase separation transition, however, occurs for 
A > 1 as given in the Fig. ([7]). In this case the transi- 
tions to phase separation and to the Mott insulator occur 
around same U c ~ 5.7. In other words we did not find 
a PS-SF phase sandwiched between 2SF and PS-MI for 
this case. 



IV. CONCLUSION 



FIG. 5: (a)Plots of gap Gl versus 1/L for different values of 
U. The gap goes to zero linearly when L — » oo for all the 
values of U considered. Here p a = 1, pb — 1/2 and A = 0.95. 
(b) Local density distribution (n°) and (n\) for p a = 1, pt — 
1/2 and A = 1.05 for two different U = 1, 4. 

(iii)Pa = l,p& = 1: 
Finally we consider double commensurate case where 
both the species of bosons undergo SF to MI phase tran- 
sition in the absence of U ab . It may be noted that for 
p a = Pb = 1, U c ~ 3.4 for U ab = 0. In Fig. ©, we 
plot the scaling of the gap LGl for A = 1.0. From this 
figure and from similar ones for A < 1 i.e., U ab < U, 
we find that the transition from 2SF to MI occurs at a 
much higher value of U = U c ~ 5.7. No SF-MI transi- 
tion observed at U ~ 3.4. Due to the collective intra and 
inter species interactions in model |T]), the species index 
become irrelevant for the phase transition. For U ab ~ U, 
the phase transition from 2SF to MI is similar to the SF- 
MI transition in single species bosons with density p = 2. 



We have studied the ground state of a two species 
Bose mixture in one dimension using the finite-size 
density-matrix renormalization group method. We 
have considered three sets of densities {p a , Pb} = 
{1/2, 1/2}, {1, 1/2}, {1, 1}. Analyzing the scaling of gap 
in the energy spectrum and the order parameter for phase 
separation we have obtained several phases: 2SF, MI, PS- 
SF and PS-MI. For U ab < U, the Mott Insulator phase is 
possible only when the total density p = p a + Pb = n, is 
an integer. The superfluid to Mott Insulator transition 
in model (TTJ) is then similar to the single species Bose- 
Hubbard model with the same total density p. The criti- 
cal on-site interaction for the 2SF-MI transition, however, 
depends on the values of A. The lower the value of A, 
the larger the value of U c . For p ^ n, Mott insulator 
phase is not found. Phase separation occurs for U ab > U 
irrespective of the value of density. For p — n and for 
all the values of A that we have considered, we found 
phase separated Mott insulator phase. In the case of 
Pa = Pb = 1/2, we observe a phase separated superfluid 
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PS-SF sandwiched between 2SF and PS-MI. However, FIG. 7: Plots of LGl (a) and Ops (b) versus U demonstrate 

for p a = pi, = 1, no PS-SF was found and the transition various phases in the case of p a = pt = 1, A = 1.05. 

is directly from 2SF to PS-MI. For p a = l,p b = 1/2 we 

found a transition from 2SF to PS-SF for A > 1 and only 

2SF phase for A < 1. It would indeed be worthwhile to 

devise experiments to test our findings. 
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